The role of dark matter interaction in galaxy clusters 
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We consider a toy model to analyze the consequences of dark matter interaction with a dark energy 
background on the overall rotation of galaxy clusters and the misalignment between their dark 
matter and baryon distributions when compared to ACDM predictions. The interaction parameters 
are found via a genetic algorithm search. The results obtained suggest that interaction is a basic 
phenomenon whose effects are detectable even in simple models of galactic dynamics. 
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I. INTRODUCTION 

Today, it is basically taken for granted that as much 
as one quarter to one third of the energy content of the 
Universe is composed by a strange and heavy kind of 
matter, the so called dark matter. An even stranger neg- 
ative pressure object, responsible for its acceleration, the 
so called dark energy, fills about two thirds. Baryons 
are responsible for less than five percent of the energy 
content of the Universe [lj . 

In general one assumes that dark energy is described by 
a cosmological constant, an assumption compatible with 
the recent WMAP data[2j. However, there is a cloud 
upon these results, since the theoretical expectation for 
a cosmological constant, if not vanishing by some un- 
known and unexpected symmetry, is 120 orders of mag- 
nitude above the observational value. Moreover, there is 
a further question mark concerning why the dark energy 
is important exactly today, or equivalently why its value 
coincides with today's energy content of the Universe. 
This question led several authors to propose a model for 
dark energy, which naturally interacts with dark matter 
1- 

It has also been realized that galaxy clusters may con- 
tain information about dark energy and dark matter in- 
teraction because the hidden sector interaction implies 
a correction to the virialization process in the cluster 
leading to quite strong constraints in the interaction pa- 
rameters [5(. 

More recently, a displacement in the angular distri- 
bution of baryon matter with respect to dark matter in 
galaxy clusters has been detected through gravitational 
lensing [6] , and it has been suspected that the alignment 
of satellite galaxies in clusters is affected by dark mat- 
ter interaction Q. In fact, hints to a contrast between 
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the ACDM results and observations in the cluster matter 
distribution may be relevant to the overall matter distri- 
bution, since Lee et al. [8] observed a departure from 
the ACDM prediction when confronted to the fact that 
the orientations of the galaxy distributions are weakly 
correlated to the dark matter distribution in the cluster. 
In Q it has been claimed that such a contrast with the 
ACDM result is consistent at a 99% confidence level. As 
this is a highly important conclusion, an explanation is 
mandatory. 

Here we follow a similar trend, showing that a displace- 
ment in the angular distribution can be traced back to 
dark matter interaction. Indeed, if dark matter interacts 
with dark energy, there will be a kind of external poten- 
tial for dark matter that does not affect baryon matter 
and the behavior of dark matter with respect to baryons 
will be unbalanced. Such a behavior has already been 
analysed in the case of galactic clusters in some recent 
papers 0, [IJ with a positive (sometimes marginal) an- 
swer for the interaction. Now we consider the effect of 
the interaction on the rotation of clusters. 

Our aim is to consider a very simple model for the 
cluster and show that an external potential mimicking 
the interaction of the dark sector is enough to lead to 
results similar to the observation. We consider the inter- 
action to be, generally speaking, such that its strength is 
proportional to the existing matter density, namely 



Pdm + 3-ffpDM = gHp , 



(1) 



where H is the Hubble constant, g is a coupling to be 
determined and p in the right-hand side is a combina- 
tion of dark matter and dark energy 0. Equation {TJ 
gives us a phenomenological effect which might account 
for the recent observations of the clusters presenting this 
anomaly. Such an interaction has been widely used to 
model the dark sector interaction, although it is not the 
only possible form. Indeed, the two-fluid interaction has 
several different possible forms and can also be originated 
from a field-theory-inspired interaction, which we do not 
pursue here. 
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II. A TOY MODEL FOR CLUSTER DYNAMICS 

There are strong evidences in cluster dynamics to en- 
courage the notion that the dark matter concentration 
has a dynamics of its own and sometimes does not fol- 
low the baryon component. In dark matter interaction, 
there is no scattering or diffusion, but only gravitational 
interaction governing the peculiar movements (mostly ro- 
tation). Individual galaxies also interact with each other 
mainly by gravitation. The classical interaction is suffi- 
cient to provide a means of separating dark matter and 
galactic baryons. However, it has been pointed out j|J 
that a purely gravitational interaction (or else, a ACDM 
scenario) does not fit the observations and the dark mat- 
ter follows a different dynamics. In such a case a non- 
gravitational force would be required. Moreover, such 
a force is highly improbable for baryons, whose interac- 
tions are too well known to allow for any new interac- 
tion. There is a quite natural new interaction for the 
dark sector in case the dark energy is not a cosmological 
constant, but a new field, presumably part of the parti- 
cle physics family, or an extended new standard model. 
The dark energy interaction with the dark matter implies 
the existence of an external potential for the dark matter 
behavior. 

One consequence of this interaction may be observed 
in astrophysical objects. It is a result from cosmologi- 
cal simulations that galaxy clusters formed from gravi- 
tational collapse are triaxial Q. After virialization, or 
sufficiently close to it, the cluster evolution may be char- 
acterized by rigid-body parameters, such as a net angular 
momentum relative to the axes. Thus, a simple model for 
the dynamics of a cluster as a whole, disregarding indi- 
vidual components, would consist of a system formed by 
two overlapping triaxial prolate ellipsoids, correspond- 
ing to the dark and baryon distributions. Their effective 
radius and angular momentum would depend on their 
density profiles, and ultimately be determined by obser- 
vation and accretion models. 

For our purposes, we go one step further on our sim- 
plifying hypotheses. We assume that the rotation is con- 
fined to only one axis. This assumption is supported by 
the fact that large mass clusters have an increased ellip- 
ticity [10]. Therefore, such a system is sufficiently well 
described by a rigid bar or a rigid set of collinear point 
masses. 

Taking an interaction such as (TTJ) as a working hy- 
pothesis, we consider this simple mechanical model for a 
cluster: a system consisting of two pairs of point masses 
mi, ?7i2, 7773 and 7774 connected by two massless poles of 
length ?b and l d , which correspond to the length of the 
semi-major axes of the baryon and dark matter distri- 
butions respectively, spinning freely around an axis fixed 
at the center of mass, interacting only through Newto- 
nian gravity (see figure [l| . Keeping the rotation axis 
fixed is merely a consequence of assuming that the dark 
and baryon components have both zero linear momen- 
tum, and not an a priori simplification. From angular 



momentum conservation, we may assume without loss of 
generality that the system is confined to the xy plane. 




FIG. 1. Toy model for cluster matter disposition. The white 
spheres (1, 2) correspond to baryon matter and the black 
spheres (3, 4) to dark matter. 

Therefore, the Lagrangian for this model can be writ- 
ten in terms of the angles between the poles and the x 
axis. Thus, 



ith 



(mi + m 2 ) l h 8 h + (m 3 + m 4 ) l d 9 d 

77717773 + 777.2TO4 mirrii + 77727773 
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ll + ll±2l h l d cos(9 h -e d ) 



D ± (6 h ,9 d ) -- 
The corresponding equations of motion are 



(2) 



(3) 
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(4a) 



(4b) 



III. 



SIMULATION PARAMETERS AND 
RESULTS 



For simplicity, we assume mi = 777 2 = mb constant and 
7773 = 7774 = m d (t). The effects of dark matter interaction 
shape the function m d (t) according to the linear model 
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m<x(t) = mo + Xt, where A = gH$. We have abandoned 
the global expansion in equation ([1} and taken the right- 
hand side to be proportional only to /9dm- We consider 
an expansion of the solution up to first order. 

We have scaled the time parameter in terms of the 
age of the clusters, such that t — 0, the starting time 
of the simulation, corresponds to their formation epoch, 
and t = 1, the final time, corresponds to today. 



A. Fixed interaction parameter 

As a preliminary approach, we have run a simulation of 
several clusters with a fixed interaction parameter g, to 
verify if there is any consequence at all to the disposition 
angles when compared to the non-interacting case where 

771d = 0. 

We have fit the parameter A = TOd(l) — too to g = 0.15 
from observation |3| so that the masses reach TOd(l) = 
0.857M T 0, with M T = m h + m d (l), with an initial 
mass given by 



TO d (0) 



TO d (l) 

1 + gH 



(5) 



where we have cast the Hubble constant in terms of the 
rescaled time. 

To ensure that the system has enough initial energy 
not to undergo gravitational collapse, the initial veloci- 
ties were chosen such that the tension on the poles would 
be zero at the initial time, therefore rendering the pole 
construction less artificial. For small A8 = (6^ — 9d), it 
may not be possible to satisfy this condition for #b(0). In 
such cases, we have adopted the same initial value com- 
puted for ^d(O), which means that the system is initially 
coupled and rotates synchronously. 

We have also cast the mass parameters in terms of ra- 
tios with respect to the baryon mass of the cluster, which 
means that we have taken mt = f, and the length pa- 
rameters as ratios with respect to the typical semi-major 
axis of the baryon component of the cluster. Consid- 
ering these factors in mass, time and length, we have 
adjusted the gravitational constant G accordingly. The 
remaining parameters used in our simulation were there- 
fore Mt = 20 /3, lb = 1 and ^d = 1-1. Using these pa- 
rameters we have integrated the system (|4"a|) . (|4bl) using a 
fourth-order Runge-Kutta method for a gaussian distri- 
bution of initial angular differences centered at 0° with 
a = 3.3° to match the angular distribution provided by 
simulations through the ACDM-based predictions 

For an initially synchronous system, as is presumably 
the case for a standard formation of a cluster from the 
accretion of both dark and baryon matter, the angu- 
lar velocities eventually decouple, and the final angular 
distance A#(l) for a starting angular difference Ada is 
roughly described by the fitting 



where A is the oscillation amplitude, ip is an unimportant 
phase and £, also in degrees, is small but non- vanishing. 
At the final time t — 1 we have observed an increased 
scatter of angle dispersions (see figure [5]). 




A6(t) = Ae st cos(wt + (p) + £ . 



(G) 
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FIG. 2. Initial angular differences and final oscillation ampli- 
tudes for a simulated cluster with parameters given in section 
Al interaction parameter g = 0.15 and 20,000 initial angle 
differences. There have been no significant occurrences above 
45°. 

For an initially Gaussian distribution, in the non- 
interacting case the final amplitude retains the initial 
shape, as can be seen by setting rhd = in the system 
(j4|). When the interaction is present, at the final time 
the distribution becomes more scattered while remaining 
spread around zero, which is in accordance with the data 



B. Genetic search for the best- fit interaction 
parameter 

For more robust comparisons, in order to find a best- 
fit value for the interaction strength g, we need an idea 
of its correlation with the angular distribution measured 
by observations. In order to achieve a figure we have run 
the program with the interaction turned on and searched 
for the best fit in comparison with the observed spread 
of the clusters. 

The search was implemented via a genetic algorithm 
whose fitness function was the % 2 computed from the 
difference between the observational data and the calcu- 
lated frequency of the angle AO in the simulated clusters. 
We have used the same initial gaussian sample of angles 
as in section MI Al We started with random values of g 
ranging from to 5. The fitness function was given by the 
X 2 between the simulated frequencies and the observed 
data from Q . 

We can use the raw data from [8| or the more narrow 
data selected in that reference. The former leads to an 
interaction constant g — 1.75 which is large when com- 
pared with other previous analyses 0, However, tak- 
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ing the more restricted set of clusters, one finds the result 
shown in figure [3J leading to an interaction strength 



g = 0.86 



(7) 



with x 2 — 1-5. The above value, although not completely 
ruling out non-interaction, fits the observed values better 
than the non-interacting case (for which the adjustment 
results in \ 2 = 1-65. Thus, there is a hint that such 
observations are compatible with a non vanishing inter- 
action in the dark sector, a conclusion which points, once 
more, towards a dark energy different from a simple cos- 
mological constant. 
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FIG. 3. Results of the best-fit search for the interaction pa- 
rameter for the low uncertainty set of observed clusters from 
for f ,000 simulated clusters, in comparison with the non- 
interacting case. 



It should be noted that, although we know that this 
very simplified Classical Mechanical model does not de- 
scribe the full dynamics of galactic clusters, the fact that 
it successfully reproduces (or better approximates) an im- 
portant feature of the observation data is a hint that dark 
matter interaction may be a very basic phenomenon with 
detectable influences even in the classical realm of inter- 
actions. 
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